0\ 

C\ '. Dynamics of Nonequilibrium Deposition 

G 

Vladimir Privman 

o 

_h ! Department of Physics, Clarkson University, Potsdam, New York 13699-5820, USA 



c3 

ABSTRACT 
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In this work we survey selected theoretical developments for models 
of deposition of extended particles, with and without surface diffusion, on 

linear and planar substrates, of interest in colloid, polymer, and certain 

i— i ' 

O ' biological systems. 
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1. Introduction 



Dynamics of important physical, chemical, and biological processes, 
e.g., [1-2], provides examples of strongly fluctuating systems in low di- 
mensions, D = 1 or 2. These processes include surface adsorption, for in- 
stance of colloid particles or proteins, possibly accompanied by diffusional 
or other relaxation (such as detachment), for which the experimentally 
relevant dimension is that of planar substrates, D = 2, or that of large 
collectors. The surface of the latter is also semi-two-dimensional owing 
to their large size as compared to the size of the deposited particles. 

For reaction-diffusion kinetics, the classical chemical studies were for 
D = 3. However, recent emphasis on heterogeneous catalysis generated 
interest in D = 2. Actually, for both deposition and reactions, some 
experimental results exist even in D = 1 (literature citations will be 
given later). Finally, kinetics of ordering and phase separation, largely 
amenable to experimental probe in D = 3 and 2, attracted much recent 
theoretical effort in D = 1, 2. 

Theoretical emphasis on low-dimensional models has been driven by 
the following interesting combination of properties. Firstly, models in 
D = 1, and sometimes in D = 2, allow derivation of analytical results. 
Secondly, it turns out that all three types of model: deposition-relaxation, 



reaction-diffusion, phase separation, are interrelated in many, but not all, 
of their properties. This observation is by no means obvious, and in fact 
it is model-dependent and can be firmly established and explored only in 
low dimensions, especially in D = 1, see, e.g., [1-2]. 

It turns out that for systems with stochastic dynamics without the 
equilibrium state, important regimes, such as the large-time asymptotic 
behavior, are frequently governed by strong fluctuations manifested in 
power-law rather than exponential time dependence, etc. However, the 
upper critical dimension above which the fluctuation behavior is described 
by the mean-field (rate-equation) approximation, is typically lower than 
in the more familiar and better studied equilibrium models. As a result, 
attention has been drawn to low dimensions where the strongly fluctuat- 
ing non-mean-field behavior can be studied. 

Low-dimensional nonequilibrium dynamical models pose several in- 
teresting challenges theoretically and numerically. While many exact, 
asymptotic, and numerical results are already available in the literature, 
as reviewed in [1-2], this field presently provides examples of properties 
(such as power-law exponents) which lack theoretical explanation even 
in ID. Numerical simulations are challenging and require large scale 
computational effort already for ID models. For more experimentally 
relevant 2D cases, where analytical results are scarce, difficulty in nu- 



merical simulations has been the "bottle neck" for understanding many 
open problems. 

The purpose of this work is to provide an introduction to the field of 
nonequilibrium surface deposition models of extended particles. No com- 
prehensive survey of the literature is attempted. Relation of deposition to 
other low-dimensional models mentioned earlier will be only referred to in 
detail in few cases. The specific models and examples selected for a more 
detailed exposition, i.e., models of deposition with diffusional relaxation, 
were biased by author's own work. 

The outline of the review is as follows. The rest of this introductory 
section is devoted to defining the specific topics of surface deposition 
to be surveyed. Section 2 describes the simplest models of random se- 
quential adsorption. Section 3 is devoted to deposition with relaxation, 
with general remarks followed by definition of the simplest, ID models of 
diffusional relaxation for which we present a more detailed description of 
various theoretical results. Multilayer deposition is also addressed in Sec- 
tion 3. More numerically-based 2D results for deposition with diffusional 
relaxation are surveyed in Section 4, along with concluding remarks. 

Surface deposition is a vast field of study. Indeed, dynamics of the 
deposition process is governed by substrate structure, substrate-particle 
interactions, particle-particle interactions, and transport mechanism of 
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particles to the surface. Furthermore, deposition processes may be ac- 
companied by particle motion on the surface and by detachment. Our 
emphasis here will be on those deposition processes where the particles 
are "large" as compared to the underlying atomic and morphological 
structure of the substrate and as compared to the range of the particle- 
particle and particle-substrate interactions. Thus, colloids, for instance, 
involve particles of submicron to several micron size. We note that 
l/xm= 10000A, whereas atomic dimensions are of order 1A, while the 
range over which particle-surface and particle-particle interactions are 
significant as compared to kT, is typically of order lOOAor less. 

Extensive theoretical study of such systems is relatively recent and 
it has been motivated by experiments where submicron-size colloid, poly- 
mer, and protein "particles" were the deposited objects; see [3-18] for a 
partial literature list, as well as other articles in this issue. It is usually 
assumed that the main mechanism by which particles "talk" to each other 
is exclusion effect due to their size. In contrast, deposition processes as- 
sociated, for instance, with crystal growth, e.g., [19], involve atomic-scale 
interactions and while the particle-particle exclusion is always an impor- 
tant factor, its interplay with other processes which affect the growth 
dynamics is quite different. 

Perhaps the simplest and the most studied model with particle ex- 



elusion is Random Sequential Adsorption (RSA). The RSA model, to be 
described in detail in Section 2, assumes that particle transport (incoming 
flux) onto the surface results in a uniform deposition attempt rate R per 
unit time and area. In the simplest formulation, one assumes that only 
monolayer deposition is allowed. This could correspond, for instance, to 
repulsive particle-particle and attractive particle-substrate forces. Within 
this monolayer deposit, each new arriving particle must either "fit in" in 
an empty area allowed by the hard-core exclusion interaction with the 
particles deposited earlier, or the deposition attempt is rejected. 

As mentioned, the basic RSA model will be described shortly, in Sec- 
tion 2. More recent work has been focused on its extensions to allow for 
particle relaxation by diffusion, Sections 3 and 4, to include detachment 
processes, and to allow multilayer formation. The latter two extensions 
will be briefly surveyed in Section 3. Many other extensions will not 
be discussed, such as for instance "softening" the hard-core interactions 
[13,20] or modifying the particle transport mechanism, etc. [21-22]. 
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2. Random Sequential Adsorption 



The irreversible Random Sequential Adsorption (RSA) process [21- 
22] models experiments of colloid and other, typically, submicron, par- 
ticle deposition [4-16] by assuming a planar 2D substrate and, in the 
simplest case, continuum (off-lattice) deposition of spherical particles. 
However, other RSA models have received attention. In 2D, noncircu- 
lar cross-section shapes as well as various lattice-deposition models were 
considered [21-22]. Several experiments on polymers [3] and attachment 
of fluorescent units on DNA molecules [18] (the latter, in fact, is usually 
accompanied by motion of these units on the DNA "substrate" and de- 
tachment) suggest consideration of the lattice-substrate RSA processes, 
in ID. RSA processes have also found applications in traffic problems 
and certain other fields and they were reviewed extensively in the litera- 
ture [21-22]. Our presentation in this section aims at defining some RSA 
models and outlining characteristic features of their dynamics. 

Figure 1 illustrates the simplest possible monolayer lattice RSA model: 
irreversible deposition of dimers on the linear lattice. An arriving dimer 
will be deposited if the underlying pair of lattice sites are both empty. 
Otherwise, it is discarded. Thus, the deposition attempt of a will suc- 
ceed. However, if the arriving particle is at b then the deposition attempt 
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will be rejected unless there is some relaxation mechanism such as de- 
tachment of dimers or monomers, or diffusional hopping. For instance, 
if c first hops to the left then later deposition of b can succeed. For d, 
the deposition is, again, not possible unless detachment and/or motion 
of monomers or whole dimers clear the appropriate landing sites marked 
by e. 

Let us consider the irreversible RSA without detachment or diffusion. 
Note that once a attaches, in Figure 1, the configuration is fully jammed 
in the interval shown. The substrate is usually assumed to be empty 
initially, at t = 0. In the course of time t, the coverage, p(t), increases 
and builds up to order 1 on the time scales of order (RV) , where R was 
defined earlier as the deposition attempt rate per unit time and "area" of 
the D-dimensional surface, while V is the particle volume. The latter is 
D-dimensional; for deposition of spheres on a planar surface, V is actually 
the cross-sections area. 

At large times the coverage approaches the jammed- state value where 
only gaps smaller than the particle size were left in the monolayer. The 
resulting state is less dense than the fully ordered "crystalline" (close- 
packed) coverage. For the D = 1 deposition shown in Figure 1 the fully 
ordered state would have p = 1. The variation of the RSA coverage is 
illustrated by the lower curve in Figure 2. 
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At early times the monolayer deposit is not dense and the deposition 
process is largely uncorrected. In this regime, mean-field like low-density 
approximation schemes are useful [23-26]. Deposition of /c-mer particles 
on the linear lattice in ID was in fact solved exactly for all times [3,27- 
28]. In D = 2, extensive numerical studies were reported [26,29-40] of the 
variation of coverage with time and large-time asymptotic behavior which 
will be discussed shortly. Some exact results for correlation properties are 
also available, in ID [27]. 

The large-time deposit has several characteristic properties that have 
attracted much theoretical interest. For lattice models, the approach to 
the jammed-state coverage is exponential [40-42]. This was shown to 
follow from the property that the final stages of deposition are in few 
sparse, well separated surviving "landing sites." Estimates of decrease in 
their density at late stages suggest that 



p(oo) - pit) ~ exp (-m D t) , (2.1) 

where £ is the lattice spacing. The coefficient in (2.1) is of order £ D /V 
if the coverage is defined as the fraction of lattice units covered, i.e., the 
dimensionless fraction of area covered, also termed the coverage fraction, 
so that coverage as density of particles per unit volume would be V~ x p. 
The detailed behavior depends of the size and shape of the depositing 
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particles as compared to the underlying lattice unit cells. 

However, for continuum off-lattice deposition, formally obtained as 
the limit £ — ► 0, the approach to the jamming coverage is power-law. 
This interesting behavior [41-42] is due to the fact that for large times the 
remaining voids accessible to particle deposition can be of sizes arbitrarily 
close to those of the depositing particles. Such voids are thus reached 
with very low probability by the depositing particles, the flux of which 
is uniformly distributed. The resulting power-law behavior depends on 
the dimensionality and particle shape. For instance, for D-dimensional 
cubes of volume V, 



p(oo) - p(t) ~ L V ^ , (2.2) 



while for spherical particles, 

p(oo) - p(t) ~ (RVt)-^ D . (2.3) 

For the linear surface, the D = 1 cubes and spheres both reduce to the 
deposition process of segments of length V. As mentioned earlier, this 
ID process is exactly solvable [27]. 

The D > 1 expressions (2.2)-(2.3), and similar relations for other 
particle shapes, etc., are actually empirical asymptotic laws which have 
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been verified, mostly for D = 2, by extensive numerical simulations [29- 
40]. The most studied 2D geometries are circles (corresponding to the 
deposition of spheres on a plane) and squares. The jamming coverages 
are [29-31,39-40] 



/Osquares(oo) ~ 0.5620 and /Ocircles(oo) ~ 0.544 to 0.550 . (2.4) 

For square particles, the crossover to continuum in the limit k — > oo and 
I — > 0, with fixed V X I D = k£ in deposition of k x k x . . . x k lattice 
squares, has been investigated in some detail [40], both analytically (in 
any D) and numerically (in 2D). 

The correlations in the large-time "jammed" state are different from 
those of the equilibrium random "gas" of particles with density near 
p(oo). In fact, the two-particle correlations in continuum deposition de- 
velop a weak singularity at contact, and correlations generally reflect the 
infinite memory (full irreversibility) of the RSA process [27,31,42]. 
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3. Deposition with Relaxation 



Monolayer deposits may "relax" (i.e., explore more configurations) 
by particle motion on the surface, by their detachment, etc. In fact, de- 
tachment has been experimentally observed in deposition of colloid parti- 
cles which were otherwise quite immobile on the surface [7]. Theoretical 
interpretation of colloid particle detachment data has proved difficult, 
however, because binding to the substrate once deposited, can be differ- 
ent for different particles, whereas the transport to the substrate, i.e., the 
flux of the arriving particles in the deposition part of the process, typi- 
cally by convective diffusion, is more uniform. Detachment also plays role 
in deposition on DNA molecules [18]. Theoretical interpretation of the 
latter data, which also involves hopping motion on DNA, was achieved 
by mean- field type modeling [43]. 

Recently, more theoretically motivated studies of the detachment re- 
laxation processes, in some instances with surface diffusion allowed as 
well, have lead to interesting model studies [44-50]. These investigations 
did not always assume detachment of the original units. For instance, in 
the ID dimer deposition shown in Figure 1, each dimer on the surface 
could detach and open up a "landing site" for future deposition. However, 
in order to allow deposition in the location represented schematically by 
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the dimer particle d, two monomers could detach (marked by e) which 
were parts of different dimers. Such models of "recombination" prior to 
detachment, of /c-mers in D = 1, were mapped onto certain spin mod- 
els and symmetry relations identified which allowed derivation of several 
exact and asymptotic results on the correlations and other properties 
[44-50]. We note that deposition and detachment combine to drive the 
dynamics into a steady state, rather than jammed state as in ordinary 
RSA. These studies have been largely limited thus far to ID models. 

We now turn to particle motion on the surface, in a monolayer de- 
posit, which was experimentally observed in deposition of proteins [17] 
and also in deposition on DNA molecules [18,43]. From now on, we focus 
on diffusional relaxation (random hopping in the lattice case). Consider 
the dimer deposition in ID; see Figure 1. The configuration in Fig- 
ure 1, after particle a is actually deposited, is jammed in the interval 
shown. Hopping of particle c one site to the left would open up a two- 
site gap to allow deposition of b. Thus, diffusional relaxation allows the 
deposition process to reach denser, in fact, ordered (close-packed) config- 
urations. For short times, when the empty area is plentiful, the effect of 
the in-surface particle motion will be small. However, for large times, the 
density will exceed that of the RSA process, as illustrated by the upper 
curve in Figure 2. 
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Further investigation of this effect is much simpler in ID than in 2D. 
Let us therefore consider the ID case first, postponing the discussion of 
2D models to the next section. Specifically, consider deposition of /c-mers 
of fixed length V. In order to allow limit k — > oo which corresponds to 
continuum deposition, we take the underlying lattice spacing £ = V/k. 
Since the deposition attempt rate R was defined per unit area (unit length 
here), it has no significant /c-dependence. However, the added diffusional 
hopping of k-mers on the ID lattice, with attempt rate H and hard- 
core or similar particle interaction, must be /c-dependent. Indeed, we 
consider each deposited /c-mer particle as randomly and independently 
attempting to move one lattice spacing to the left or to the right with 
rate H/2 per unit time. Of course, particles cannot run over each other 
so some sort of hard-core interaction must be assumed, i.e., in a dense 
state most hopping attempts will fail. However, if left alone, each particle 
would move diffusively for large time scales. In order to have the resulting 
diffusion constant T> finite in the continuum limit k — > oo, we put 



H oc V/£ 2 = Vk 2 /V 2 . (3.1) 

which is only valid in ID. 

Each successful hopping of a particle results in motion of one empty 
lattice site (see particle c in Figure 1). It is useful to reconsider the dy- 
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namics of particle hopping in terms of the dynamics of this rearrangement 
of empty area fragments [51-53]. Indeed, if several of these empty sites are 
combined to form large enough voids, deposition attempts can succeed 
in regions of particle density which would be "frozen" or "jammed" in 
ordinary RSA. In terms of these new "particles" which are empty lattice 
sites of the deposition problem, the process is in fact that of reaction- 
diffusion. Indeed, k reactants (empty sites) must be brought together by 
diffusional hopping in order to have finite probability of their annihila- 
tion, i.e., disappearance of a group of consecutive nearest-neighbor empty 
sites due to successful deposition. Of course, the /c-group can also be bro- 
ken apart due to diffusion. Therefore, the /c-reactant annihilation is not 
instantaneous in the reaction nomenclature. Such /c-particle reactions are 
of interest on their own [54-59]. 

The simplest mean-field rate equation for annihilation of k reactants 
describes the time dependence of the coverage, p(t), in terms of the reac- 
tant density 1 — p, 



| = r(i - „)* , 0.2) 



where T is the effective rate constant. Note that we assume that the close- 
packing dimensional coverage is 1 in ID. There are two problems with 
this approximation. Firstly, it turns out that for k = 2 the mean-field 

- 15- 



approach breaks down. Diffusive-fluctuation arguments for non-mean- 
field behavior have been advanced for reactions [54,56,60-61]. In ID, 
several exact calculations support this conclusion [62-68] . The asymptotic 
large-time behavior turns out to be 



l-p~l/Vt (k = 2,D = l) , (3.3) 

rather than the mean-field prediction ~ 1/t. The coefficient in (3.3) is 
expected to be universal, when expressed in an appropriate dimensionless 
form by introducing single-reactant diffusion constant. 

The power law (3.3) was confirmed by extensive numerical simula- 
tions of dimer deposition [69] and by exact solution for one particular 
value of H [70] for a model with dimer dissociation. The latter work 
also yielded some exact results for correlations. Specifically, while the 
connected particle-particle correlations spread diffusively in space, their 
decay it time is nondiffusive; see [70] for details. Series expansion stud- 
ies of models of dimer deposition with diffusional hopping of the whole 
dimers or their "dissociation" into hopping monomers, has confirmed the 
expected asymptotic behavior and also provided estimates of the coverage 
as a function of time [71]. 

The case k = 3 is marginal with the mean-field power law modified 
by logarithmic terms. The latter were not observed in Monte Carlo stud- 
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ies of deposition [52]. However, extensive results are available directly 
for three-body reactions [56-59], including verification of the logarithmic 
corrections to the mean-field behavior [57-59]. 

The second problem with the mean-field rate equation was identified 
in the continuum limit of off-lattice deposition, i.e., for k — > oo. In- 
deed, the mean-field approach is essentially the fast diffusion approxima- 
tion assuming that diffusional relaxation is efficient enough to equilibrate 
nonuniform density profile fluctuations on the time scales fast as com- 
pared to the time scales of the deposition events. Thus, the mean-field 
results are formulated in terms of the uniform properties, such as the den- 
sity. It turns out, however, that the simplest, /c th -power of the reactant 
density form (3.2) is only appropriate for times t » e k ~ 1 /(RV). 

This conclusion was reached [51] by assuming the fast-diffusion, ran- 
domized (equilibrium) hard-core reactant system form of the inter-reactant 
distribution function in ID (essentially, an assumption on the form of 
certain correlations). This approach, not detailed here, allows Ginzburg- 
criterion-like estimation of the limits of validity of the mean-field results 
and it correctly suggests mean-field validity for k = 4, 5, . . ., with loga- 
rithmic corrections for k = 3 and complete breakdown of the mean-field 
assumptions for k = 2. However, this detailed analysis yields the modified 
mean-field relation 
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dp jRV(l - p)K 

dt (i -p + fc-ip) ^ V ' {AA > 

where 7 is some effective dimensionless rate constant. This new expres- 
sion applies uniformly as k — > 00. Thus, the continuum deposition is also 
asymptotically mean-field, with the essentially-singular "rate equation" 



^ = 7(1 - p) exp[-p/(l - p)] (k = 00, D = 1) . (3.5) 

The approach to the full, saturation coverage for large times is extremely 
slow, 



Similar predictions for /c-particle reactions can be found in [55]. 

When particles are allowed to attach also on top of each other, with 
possibly some rearrangement processes allowed as well, multilayer de- 
posits will be formed. It is important to note that the large-layer struc- 
ture of the deposit and fluctuation properties of the growing surface will 
be determined by the transport mechanism of particles to the surface and 
by the allowed relaxations (rearrangements). Indeed, these two charac- 
teristics determine the screening properties of the multilayer formation 
process which in turn shape the deposit morphology, which can range 
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from fractal to dense, and the roughening of the growing deposit sur- 
face. There is a large body of research studying such growth, with recent 
emphasis on the growing surface fluctuation properties. 

However, the feature characteristic of the RSA process, i.e., the ex- 
clusion due to particle size, plays no role in determining the universal, 
large-scale properties of "thick" deposits and their surfaces. Indeed, the 
RSA-like jamming will be only important for detailed morphology of the 
first few layers in a multilayer deposit. However, it turns out that RSA- 
like approaches (with relaxation) can be useful in modeling granular com- 
paction [72]. 

In view of the above remarks, multilayer deposition models involving 
jamming effects were relatively less studied. They can be divided into two 
groups. Firstly, structure of the deposit in the first few layers is of interest 
[73-75] since they retain "memory" of the surface. Variation of density 
and other correlation properties away from the wall has structure on the 
length scales of particle size. These typically oscillatory features decay 
away with the distance from the wall. Numerical Monte Carlo simulation 
aspects of continuum multilayer deposition (ballistic deposition of 3D 
balls) were reviewed in [75]. Secondly, few-layer deposition processes 
have been of interest in some experimental systems. Mean-field theories 
of multilayer deposition with particle size and interactions accounted for 
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were formulated [76] and used to fit such data [12,14-16], 
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4. Two-Dimensional Deposition with Diffusional Relaxation 



We now turn to the 2D case of deposition of extended objects on pla- 
nar substrates, accompanied by diffusional relaxation (assuming mono- 
layer deposits). We note that the available theoretical results are limited 
to few studies [38,77-79]. They indicate a rich pattern of new effects 
as compared to ID. In fact, there exists extensive literature, e.g., [81] 
on deposition with diffusional relaxation in other models, in particular 
those where the jamming effect is not present or plays no significant role. 
These include, e.g., deposition of "monomer" particles which align with 
the underlying lattice without jamming, as well as models where many 
layers are formed (mentioned in the preceding section). 

The 2D deposition with relaxation of extended objects is of interest 
in certain experimental systems where the depositing objects are proteins 
[17]. Here we focus on the combined effect of jamming and diffusion, and 
we emphasize dynamics at large times. For early stages of the deposi- 
tion process, low-density approximation schemes can be used. One such 
application was reported in [38] for continuum deposition of circles on a 
plane. 

In order to identify features new to 2D, let us consider deposition 
of 2 x 2 squares on the square lattice. The particles are exactly aligned 
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with the 2x2 lattice sites as shown in Figure 3. Furthermore, we assume 
that the diffusional hopping is along the lattice directions ±x and ±y, 
one lattice spacing at a time. In this model dense configurations involve 
domains of four phases as shown in Figure 3. As a result, immobile frag- 
ments of empty area can exist. Each such single-site vacancy (Figure 3) 
serves as a meeting point of four domain walls. By "immobile" we mean 
that the vacancy cannot move due to local motion of the surrounding 
particles. For it to move, a larger empty-area fragment must first arrive, 
along one of the domain walls. One such larger empty void is shown in 
Figure 3. Note that it serves as a kink in the domain wall. 

Existence of immobile vacancies suggests possible "frozen," glassy 
behavior with extremely slow relaxation, at least locally. In fact, the 
full characterization of the dynamics of this model requires further study. 
The first numerical results [77] do provide some answers which will be 
reviewed shortly. We first consider a simpler model depicted in Figure 4. 
In this model [78-79] the extended particles are squares of size V2 x y/2. 
They are rotated 45° with respect to the underlying square lattice. Their 
diffusion, however, is along the vertical and horizontal lattice axes, by 
hopping one lattice spacing at a time. The equilibrium variant of this 
model (without deposition, with fixed particle density) is the well-studied 
hard-square model [82] which, at large densities, phase separates into two 
distinct phases. These two phases also play role in the late stages of RSA 
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with diffusion. Indeed, at large densities the empty area is stored in 
domain walls separating ordered regions. One such domain wall is shown 
in Figure 4. Snapshots of actual Monte Carlo simulation results can be 
found in [78-79]. 

Figure 4 illustrates the process of ordering which essentially amounts 
to shortening of domain walls. In Figure 4, the domain wall gets shorter 
after the shaded particles diffusively rearrange to open up a deposition 
slot which can be covered by an arriving particle. Numerical simulations 
[78-79] find behavior reminiscent of the low-temperature equilibrium or- 
dering processes [83-85] driven by diffusive evolution of the domain-wall 
structure. For instance, the remaining uncovered area vanishes according 
to 



1 - p(t) ~ -^ . (4.F 



This quantity, however, also measures the length of domain walls in the 
system (at large times). Thus, disregarding finite-size effects and as- 
suming that the domain walls are not too convoluted (as confirmed by 
numerical simulations), we conclude that the power law (4.1) corresponds 
to typical domain sizes growing as ~ y/t, reminiscent of the equilibrium 
ordering processes of systems with nonconserved order parameter dynam- 
ics [83-85]. 
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We now turn again to the 2x2 model of Figure 3. The equilibrium 
variant of this model corresponds to hard-squares with both nearest and 
next-nearest neighbor exclusion [82,86-87]. It has been studied in lesser 
detail than the two-phase hard-square model described in the preceding 
paragraphs. In fact, the equilibrium phase transition has not been fully 
classified (while it was Ising for the simpler model). The ordering at 
low temperatures and high densities was studied in [86]. However, many 
features noted, for instance large entropy of the ordered arrangements, 
require further study. The dynamical variant (RSA with diffusion) of this 
model was studied numerically in [77]. The structure of the single-site 
frozen vacancies and associated network of domain walls turns out to 
be boundary-condition sensitive. For periodic boundary conditions the 
density "freezes" at values 1 — p ~ L~ l , where L is the linear system size. 

Preliminary indications were found [77] that the domain size and 
shape distributions in such a frozen state are nontrivial. Extrapolation 
L — ► oo indicates that the power law behavior similar to (4.1) is nondiffu- 
sive: the exponent 1/2 is replaced by ~ 0.57. However, the density of the 
smallest mobile vacancies, i.e., dimer kinks in domain walls, one of which 
is illustrated in Figure 3, does decrease diffusively. Further studies are 
needed to fully clarify the ordering process associated with the approach 
to the full coverage as t — > oo and L — > oo in this model. 
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Even more complicated behaviors are possible when the deposit- 
ing objects are not symmetric and can have several orientations as they 
reach the substrate. In addition to translational diffusion (hopping), one 
has to consider possible rotational motion. The square-lattice deposition 
of dimers, with hopping processes including one lattice spacing motion 
along the dimer axis and 90° rotations about a constituent monomer, was 
studied in [80]. The dimers were allowed to deposit vertically and hori- 
zontally. In this case [80] the full close-packed coverage is not achieved 
at all because the frozen vacancy sites can be embedded in, and move by 
diffusion in, extended structures of different "topologies." These struc- 
tures are probably less efficiently "demolished" by the motion of mobile 
vacancies than the elimination of localized frozen vacancies in the model 
of Figure 3. 

In summary, we reviewed the deposition processes involving extended 
objects, with jamming and its interplay with diffusional relaxation yield- 
ing interesting new dynamics of approach to the large-time state. While 
significant progress has been achieved in ID, the 2D systems require fur- 
ther investigations. Mean-field and low-density approximations can be 
used in many instances for large enough dimensions, for short times, and 
for particle sizes larger than few lattice units. Added diffusion allows for- 
mation of denser deposits and leads to power-law large-time tails which, 
in ID, were related to diffusion-limited reactions, while in 2D, associ- 
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ated with evolution of domain- wall network and defects, reminiscent of 
equilibrium ordering processes. 
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Figure Captions 



Figure 1: Deposition of dimers on the ID lattice. Once the arriving 
dimer a attaches, the configuration shown will be fully jammed in the 
interval displayed. Further deposition can only proceed if dimer or 
monomer diffusion (hopping) and/or detachment are allowed. Letter 
labels 6, c, d are referred to in the text. 



Figure 2: Schematic variation of the coverage fraction p(t) with time 
for lattice deposition without (lower curve) and with (upper curve) 
diffusional or other relaxation. The "ordered" density corresponds 
to close packing. Note that the short-time behavior deviates from 
linear at times of order 1/(RV). (Quantities R, V are defined in the 
text.) 



Figure 3: Fragment of a deposit configuration in the deposition of 2 x 2 
squares. Illustrated are one single-site frozen vacancy at which four 
domain walls converge (indicated by heavy lines), as well as one 
dimer vacancy which causes a kink in one of the domain walls. 
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Figure 4: Illustration of deposition of V2 x y/2 particles on the square 
lattice. Diffusional motion during time interval from t\ to £2 can 
rearrange the empty area "stored" in the domain wall to open up 
a new landing site for deposition. This is illustrated by the shaded 
particles. 
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